Biotechnology for Biofuels

Glucoamylase is essential for industrial starch saccharification, but the limited thermostability and near-neutral pH tolerance of fungal glucoamylases necessitate cooling and acidification of liquefied starch. Here, we developed an artificial intelligence-guided strategy to simultaneously improve the thermostability, pH tolerance, and catalytic activity of glucoamylase from Penicillium oxalicum (PoGA). Two property-specific machine-learning models, CASPE-T and CASPE-A, identified substitutions associated with thermostability and pH tolerance, respectively. Experimental screening identified beneficial substitutions in 11 of 21 CASPE-T and 12 of 22 CASPE-A candidates. Folding-energy-guided recombination integrated the two traits while maintaining structural compatibility. The optimal variant, PoGA T513E/Q305N, exhibited 2.21-fold higher specific activity than the wild type, with half-life extended from 22.3 to 57.9 min at 60 degrees C and from 16.6 to 64.7 min at pH 8.0. Molecular dynamics simulations attributed these improvements to reinforcement of high-occupancy hydrogen-bonding networks, suppression of conformational fluctuations in the linker and carbohydrate-binding module, enhanced long-range dynamic coordination, and preservation of a compact catalytic architecture. At 60 degrees C and pH 6.5 without acidification, PoGA T513E/Q305N produced 219.9 g/L glucose and achieved 89.1% starch conversion, 31.4% higher than the wild type. This work provides an efficient framework for multi-objective enzyme engineering and sustainable starch biorefining.

Synthetic allopolyploidization is a powerful strategy for generating strains with novel industrially relevant traits. During the design of select, hybridize, evolved and learn cycle, the mitochondrial inheritance is considered neutral and most studies overcome their phenotypic impact. Here, we generated 168 synthetic allotetraploid strains representing all pairwise parental combinations of eight wild Saccharomyces species, and systematically controlled mitochondrial inheritance in these crosses. These allopolyploids were phenotyped across 31 environmental conditions relevant to biotechnology. In ~50% of all condition-strain combinations conditions, allotetraploids exhibited mid-parent heterosis. In approximately 25% of conditions, we found crosses with hybrid vigor, particularly in combinations involving Saccharomyces arboricola x Saccharomyces mikatae . However, allotetraploids were highly unstable and rapidly sporulated. Our analyses revealed that environmental conditions accounted for the largest effects, follow by its interaction with mitotype. However, mitotype alone showed small effects but significantly affected 42% of the tested conditions, underscoring the importance of mitochondrial inheritance as a key determinant of allopolyploid performance. These findings demonstrate that mitochondrial inheritance is a key, environment-dependent determinant of allopolyploid fitness and should be considered alongside nuclear genome composition in the rational design of industrial allopolyploid strains

Ugba, a traditional alkaline-fermented condiment produced from African oil bean seeds (Pentaclethra macrophylla Benth.), remains poorly characterised at genome resolution. This study characterised a single spontaneous ugba fermentation using genome-resolved shotgun metagenomics across Early (0-24 h), Mid (48-96 h), and Late (120-144 h) composite stages. Thirty metagenome-assembled genomes (MAGs) were recovered and grouped according to their temporal abundance patterns. The Early-stage retained MAG population was dominated by Bacteroidota-associated MAGs, which exhibited the highest carbohydrate-active enzyme density and broad glycoside hydrolase repertoires, consistent with substantial genome-encoded potential for utilisation of the plant-derived seed matrix. The Mid stage showed marked restructuring of MAG composition, while Core and Late-associated MAGs increased towards fermentation maturation. Functional differentiation was also evident among Corynebacterium MAGs: C. nuruki encoded a complete urease system, whereas C. phoceense_A encoded glutamate dehydrogenase and nitrate-reduction genes. These stage-associated changes coincided with progressive alkalinisation from pH 7.15 to 8.01. Together, the results show that temporal restructuring of reconstructed microbial populations was accompanied by differentiation in genome-encoded functional potential. This study provides a genome-resolved characterisation of microbial composition and functional potential during spontaneous ugba fermentation and expands current understanding of the microbial organisation of this traditional African fermented food.

Microalgae are promising biological platforms for sustainable energy production because they can convert photosynthetic activity into electricity in biophotovoltaic (BPV) systems. Identifying strains that grow well at elevated temperatures and under reduced irradiance is an important step toward BPV applications in seasonally variable environments. Microalgae were isolated from Seondong Lake (Seoul, Republic of Korea) during three sampling periods and screened for growth at different temperatures. Among 45 isolates, KHUC-11 showed the highest weighted growth index based on growth rates at 25, 30, and 35 under dark conditions and was selected for further characterization. Phylogenetic reconstruction of the 45S rRNA gene cluster, internal transcribed spacer 2 (ITS2) secondary-structure comparison, and morphological characterization supported the description of KHUC-11 as Desmodesmus ecotolerans sp. nov. D. ecotolerans produced more biomass than Chlamydomonas reinhardtii CC-125 under both moderate- and low-light conditions. Under low light, its final dry biomass concentrations were 1.81 and 1.58 times those of CC-125 under mixotrophic and photoautotrophic conditions, respectively. Despite its smaller cell size, D. ecotolerans showed higher photosynthetic performance than CC-125, and screen-printed electrode (SPE)-based electrochemical assays demonstrated enhanced anodic current generation and CV-derived estimated power output. Taken together, these findings identify D. ecotolerans as a novel climate-resilient Desmodesmus species with bioelectrochemical properties that support further evaluation for BPV applications.

The valorisation of biological polymers requires microbial systems that can both access recalcitrant substrates and convert the resulting carbon into useful products. Although microbial genome and metagenome resources provide an expanding reservoir of candidate depolymerizing and modifying enzymes, most discovery workflows remain disconnected from enzyme optimisation and host adaptation. Here we present a coupled sequence-based enzyme discovery, enzyme evolution and synthetic yeast chassis adaptation strategy for microbial biopolymer valorisation. Focusing on laccases for the depolymerisation of lignin as a proof of concept, we combine sequence data mining for enzyme discovery, modular yeast surface display for functional screening, directed evolution for enzyme optimisation and synthetic yeast genome diversification for chassis improvement. In our study, surface display enabled functional benchmarking and recovery of improved laccase variants and synthetic-genome-enabled diversification provided a route to explore host configurations that influence display and enzyme performance. By integrating enzyme-level and chassis-level optimisation, this framework addresses a central bottleneck in converting microbial biodiversity by computational sequence repository mining into deployable biomanufacturing systems. Our results establish laccases as tractable entry points for oxidative biopolymer conversion and provide a generalizable platform for engineering yeast systems for sustainable carbon valorisation.

Genome transplantation is a key technology for synthetic genomics, enabling entire genomes to be transferred into recipient cells. Despite its importance, genome transplantation remains confined to a small number of wall-less Mollicute species and is often characterized by low efficiencies, limiting the development and testing of synthetic genomes. Because genome transplantation relies on the same polyethylene glycol mediated DNA delivery process used for plasmid transformation, improving DNA uptake is an important step toward more efficient transplantation systems. Here, we systematically optimized polyethylene glycol-mediated transformation in Mesoplasma florum, a fast-growing, non-pathogenic Mollicute with a reduced genome that represents an attractive chassis for synthetic genomics. We evaluated 20 parameters spanning cell physiology, DNA preparation, membrane conditioning, and recovery conditions. DNA topology emerged as the strongest determinant of transformation efficiency, with highly compacted DNA preparations producing up to three orders of magnitude more transformants than conventional plasmid preparations. Growth phase and polyethylene glycol concentration also strongly influenced transformation outcomes. To identify potential recipients for future genome transplantation experiments, we further screened ten strains belonging to the Mesoplasma lineage. Transformation efficiencies varied widely among strains, and Mesoplasma entomophilum W17 emerged as a particularly promising candidate. Moreover, W17 supported replication of plasmids carrying the M. florum oriC, suggesting compatibility between the replication systems of the two species. Together, these results establish an improved transformation workflow for M. florum, identify DNA topology as a major determinant of DNA uptake, and reveal a new candidate recipient for future genome transplantation studies. These advances provide practical tools and biological insights for synthetic genome engineering in Mollicutes.

Electromicrobial production (EMP) technologies aim to combine renewable electricity, CO2, and engineered microbes to make energy-dense molecules at efficiencies exceeding photosynthesis. CO2 can be electrochemically reduced to formate, which is far easier to handle at the bench than H2 or an electrode, but formate carries only two electrons per carbon against the six in a biofuel. The remaining electrons must come from oxidizing additional formate, from H2 oxidation, or from extracellular electron uptake (EEU), and no rigorous comparison of these options coupled to the choice of carbon assimilation pathway currently exists. We calculate upper-limit efficiencies for butanol production by six carbon assimilation pathways, each paired with all three electron delivery mechanisms, using electrochemical parameters drawn from a survey of the recent literature. Electrical to butanol energy conversion efficiencies range from 35.5 to 51.7%, corresponding to solar-to-fuel efficiencies of 11.7 to 17%, so even the least efficient route exceeds the 8% theoretical ceiling of algal photosynthesis. The serine variant of the reductive glycine pathway reaches an electrical energy conversion efficiency of when using H2 oxidation, within 1.9 points of the most efficient pathway, and is the only high-efficiency option that tolerates O2. This makes an EMP system that combines electron delivery by formate coupled with the serine variant of reductive glycine pathway highly attractive, as it presents few barriers to rapid, iterative engineering in the lab, and a high theoretical ceiling. Drawing both carbon and electrons from formate costs 6.2 points against H2 at a state-of-the-art whole-cell voltage (2.2 V). However, this small penalty is amplified three-fold by any rise in the CO2-to-formate cell voltage, and reaches 13.5 points at the highest whole-cell voltages reported for scaled-up CO2-to-formate electrolyzers, where formate-only operation falls to 11.2 electrical-to-fuel and 3.7% solar-to-fuel efficiency, below the ceiling of photosynthesis, against 24.7 and 8.1% for H2 (only just above algal photosynthesis). Our choice between a formate-only system and one coupled to H2 oxidation or EEU therefore depends on our belief about the trajectory of CO2 reduction technology. If whole-cell voltages continue to fall at the rate of the past decade, formate alone is the right target, and the simplicity of its workflow is bought at low cost. However, if that improvement plateaus, the electron delivery mechanism must be swappable, and a system should be designed from the outset so that it can be. At the US Department of Energy SunShot target of 2 cents per kilowatt hour, the electricity to make a US gallon of butanol costs $1.40 for a formate-only system at the state of the art, rising to $5.45 at the highest scaled-up electrolyzer voltage reported, against $1.23 and $2.47 for formate and H2 system.

The unprecedented rise in greenhouse gases such as carbon dioxide (CO2), and their detrimental effects on the atmosphere have intensified the current climate emergency. This challenge has inspired the development of novel carbon capture and utilisation (CCU) technologies, specifically microbial electrosynthesis (MES), for the bioelectrochemical fixation of CO2 into commodity chemical compounds (CCCs) using microbial biocatalysts. In this study, we systematically enriched and maintained robust mixed microbial communities of electroactive bacteria (EAB) from wastewater treatment sludge to serve as biocatalysts for MES. The MES performance of these biocatalysts under different operational conditions in bioelectrochemical systems (BES) reactors were investigated. BES reactors with enriched biocatalysts incubated and maintained under controlled conditions produced higher CCCs yields than those developed at ambient temperatures, while strict anaerobic operations further improved the MES efficiency compared to aerobic conditions. At an applied potential of -1000 mV vs Ag/AgCl, CO2 conversion reached up to 88.11%, significantly outperforming the lower potential (-600 mV) operation. These findings underscore the importance of systematic development of robust and stable biocatalysts, as well as maintaining strict anaerobic conditions, for achieving efficient MES performance, thereby demonstrating the potential of our approach for scalable industrial CCU applications.

Engineering bioenergy crops for enhanced biomass digestibility and high-value co-products is a promising approach to sustainable biofuel production. However, the impact of these engineered crops on the belowground microbiome remains poorly understood. In this study, we characterized the composition and functional responses of the belowground microbiome associated with transgenic sorghum lines: one with reduced cell wall crosslinking (AT10) and another with elevated protocatechuic acid (PCA) production (Qsub). The lines were grown under field conditions with four irrigation levels, and samples were collected twice during the growing season. We found that the two transgenic types drove a subtle but significant divergence in bacterial communities, both compositionally and functionally. Community-wide functional profiling showed that AT10 rhizospheres were consistently depleted in genes associated with mobile genetic elements, while Qsub rhizospheres were enriched in TonB-dependent transport and vanillate monooxygenase, an enzyme that funnels lignin-derived monomers toward PCA degradation, alongside reduced pectin- and cellulose-degradation genes. Genome-resolved analysis linked these community-level shifts to specific taxa: in Qsub, the enriched MAGs were taxonomically diverse, ranging from genome-reduced Patescibacteria with an episymbiotic lifestyle to taxa enriched in lignin-associated auxiliary activity enzymes. In AT10, the depleted Azospirillaceae MAGs carried unusually high mobile genetic element content, consistent with the community-wide pattern. Collectively, these findings offer new insights into how genetic modifications aimed at improving bioenergy traits shapes the composition and function of the rhizosphere microbiome, with potential downstream consequences for soil carbon cycling.

Extremophilic fungi are a relatively unexplored source of stress-resilient microbes for sustainable agriculture. Here, we describe Talaromyces pseudofuniculosus FAR, isolated from the highly acidic site of Piccolo Inferno in southern Italy and identified by ITS and {beta}-tubulin (BenA) sequencing. Despite its origin from an extreme environment, T. pseudofuniculosus FRA remained fully functional under non-stress conditions and displayed tolerance to salinity and heavy metals while maintaining key plant growth-promoting traits, including phosphate solubilization, auxin-related compound production and hydrolytic activities. In soybean (Glycine max L.), the live fungus promoted plant growth under both non-stress conditions and agriculturally relevant salinity and copper stress, and successfully colonized the roots. A cell-free extracellular extract (exCFS) promoted early seedling development in vitro, but its effects were not sustained in pot experiments. GC-MS analysis provided a qualitative characterization of the extracellular metabolite profile and revealed condition-dependent differences between control, NaCl and CuSO4; growth conditions, indicating environment-dependent modulation of fungal metabolism. Overall, these findings highlight the physiological and metabolic plasticity of T. pseudofuniculosus FRA and support its potential as a fungal biostimulant for saline- and metal-impacted agricultural soils.

Anaerobic gut fungi (AGF) are an integral part of the herbivore gut microbiome, primary decomposers of lignocellulose, and play a crucial role in establishing foundational links of syntrophic associations. The syntrophy with methanogens appears beneficial for mitigating feedback inhibition and results in the augmentation of AGF activity. However, the thermodynamic advantages of the presence of methanogens to AGF remain poorly understood. This study explores the preferred pathways adopted by Cyllamyces sp. and Methanobrevibacter olleyae, using 13C-glucose as a substrate and spiked with malate and pyruvate in different experimental set-ups. The results indicated reduced biomass formation by Cyllamyces sp. in monoxenic culture with Methanobrevibacter olleyae compared to axenic AGF culture supplemented with glucose. Additionally, metabolic analyses revealed the pyruvate:ferredoxin oxidoreductase (PFO) as a key contributor to pyruvate metabolism. Malate supplementation notably enhanced the growth of Cyllamyces sp. and M. olleyae, whereas pyruvate supplementation enhanced the growth of Cyllamyces sp. only in axenic culture. This investigation provides insights into the preferred metabolic pathways operating in the hydrogenosomes of axenic anaerobic gut fungi (AGF) cultures and monoxenic AGF-methanogen co-cultures.

L-Threonine is an essential amino acid widely used in the feed, food, and pharmaceutical industries, and microbial fermentation is now the dominant route for its production. Further increases in producer titers depend on random mutagenesis coupled with high-throughput screening (HTS), but progress has been limited by the lack of a biosensor with a broad operating range, high sensitivity, and target specificity. Here, we show that L-threonine induces the tryptophan-repressed trp promoter of Escherichia coli, a response that, to our knowledge, has not been described previously. Deletion analysis showed that both the trpO operator and the trpL attenuation element contribute to induction, and that the lowest basal expression occurs only when both elements are present. By optimizing the ribosome-binding site, the promoter -10 region, and TrpR expression, we constructed an L-threonine biosensor on a pACYC backbone. The sensor responded linearly to L-threonine from 0.1 to 30 g/L (R2 = 0.9895), with an induction ratio of up to 51.7-fold. Among the 20 proteinogenic amino acids tested, only L-threonine produced a detectable response. The biosensor distinguished strains with different L-threonine production capacities and guided the directed evolution of homoserine kinase (ThrB) at a conserved motif, yielding variants with brighter reporter fluorescence and increased L-threonine production. This work provides a high-performance biosensor based on native E. coli regulatory elements for screening L-threonine overproducers.

High temperatures during grain filling stage are becoming increasingly frequent, compromising both grain and eating quality and thereby driving demand for heat-resilient cultivars. Such conditions are known to reduce the expression of granule-bound starch synthase I (GBSSI) and starch branching enzyme IIb (BEIIb), which are involved in starch biosynthesis, resulting in a decrease in amylose content and an increase in long-chain amylopectin. Thus, the present study introduced functional mutation in starch synthase IIa (SSIIa) that increases the proportion of short amylopectin chains to genetically compensate for the high-temperature-induced increase in amylopectin long chain. Rice lines carrying the ss2a mutation were grown at two locations with cooler (Akita) and warmer (Okayama) temperatures. Their grain traits, starch structure, and eating quality were compared. The ss2a mutant lines showed an increased proportion of short amylopectin chains as well as an increased apparent amylose content. Furthermore, these alterations in starch structure varied with the grain-filling temperature of the cultivation sites, ultimately affected eating quality. These results suggest that enriching short amylopectin chain via the ss2a mutation can counteract the increase in long amylopectin chain caused by high temperatures during grain filling, thereby maintaining a desirable starch structure and eating quality.

Bacterial genomes are dynamic and harbour accessory regions that can be gained, lost, or rearranged in response to selective pressures. Members of the genus Komagataeibacter are used for bacterial cellulose production and researched as platforms for engineered living materials. However, their genomes are rich in mobile genetic elements (MGEs), including insertion sequences (ISs), genomic islands and prophage remnants, which drive genome rearrangements under agitated cultivation conditions. Whether this genome architecture contributes to genomic changes during routine strain maintenance and distribution remains poorly understood. Comparative genome sequencing of independent Komagataeibacter rhaeticus iGEM stocks revealed recurrent loss of a 210 kbp chromosomal region flanked by homologous IS5 sequences. This deletion was identified in a stock received as an agar stab and reproduced from an intact-genome population subjected to simulated agar-stab handling. Phylogenomic analysis of the affected locus revealed gene-tree/species-tree discordance, heterogeneous GC content and multiple ISs, supporting its characterization as a dynamic accessory region. Populations carrying the larger genome reduction showed improved growth and cellulose production relative to the intact genome reference, whereas populations carrying a smaller nested deletion showed reduced performance. These findings demonstrate that routine strain handling is not necessarily genomically neutral and can expose structural genomic variabilities that generate population-level genomic heterogeneity with phenotypic consequences. For bacteria with MGE-rich genomes, handling history may represent an underappreciated source of genomic and phenotypic variation between laboratories. More broadly, identifying unstable accessory regions may highlight genomic content that is dispensable under laboratory conditions, providing candidates for genome streamlining and strain improvement.

Acinetobacter sp. strain Tol 5 is an attractive microorganism because of its ability to utilize various aromatic compounds, including toluene. However, because Tol 5 is a non-model microorganism, the utilization of carbon sources in central carbon metabolism remains unclear. In this study, the activity of the central carbon metabolism was characterized using 13C-metabolic flux analysis (13C-MFA). Tol 5 was cultured in synthetic medium using ethanol or acetate as the sole carbon source. In addition, because toluene is degraded into pyruvate and acetyl-CoA via the toluene dioxygenase pathway, this situation was mimicked by culturing with ethanol and lactate as co-substrates. The 13C-MFA confirmed flux through the glyoxylate shunt, which bypasses part of the TCA cycle when either ethanol or acetate was used as the sole carbon source. In ethanol, the flux ratio of isocitrate dehydrogenase at the isocitrate branch point decreased markedly. Unlike acetate, the conversion of ethanol to acetyl-CoA is accompanied by NADH production. Therefore, the observed change in this flux ratio was influenced by the excess NADH generated. Under ethanol-lactate conditions, Tol 5 simultaneously utilized ethanol and lactate. The ethanol consumption rate decreased to less than half of that observed when ethanol was the sole carbon source. Most of the lactate uptake was converted to pyruvate and metabolized via the TCA cycle. The glyoxylate shunt flux was lower than that in ethanol. These metabolic characteristics provide valuable insights into the use of Tol 5 as a next-generation microbial cell factory.

Azetidine-2-carboxylic acid (Aze) is a toxic, non-proteogenic amino acid produced in restricted plants. It can enter the human food chain through direct or indirect dietary exposure and has been linked to neurological disorders due to its ability to misincorporate for its structural analog proline (Pro) during translation and lead to protein misfolding. Although several microorganisms have evolved mechanisms to detoxify Aze, these pathways have not previously been exploited to engineer Aze tolerance in plants. Here, we heterologously expressed a bacterial Aze hydrolase (HAS) in Arabidopsis thaliana to establish a metabolic strategy for Aze detoxification in plants. Transgenic plants exhibited markedly enhanced tolerance to Aze, as demonstrated by improved root growth, cotyledon development, and survival compared with wild-type plants under Aze stress. We further showed that HAS localized to the cytosol, where it converted Aze into the non-toxic metabolite 2-hydroxy-4-aminobutyrate (HAB). HAB formation was confirmed both in planta and in vitro, demonstrating that HAS remained catalytically active in plant cells and effectively detoxified Aze. Enzyme kinetics revealed that HAS has a high affinity for Aze and is stereospecific for L-Aze, with no activity toward D-Aze. Multi-omics analysis including transcriptomics and misincorporation proteomics show that the transgenic lines misincorporate Aze significantly less compared to wild-type leading to enhanced Aze tolerance. Our findings establish bacterial HAS as an effective metabolic detoxification system that confers enhanced Aze tolerance in plants. This study provides a promising strategy for engineering plant resistance to toxic non-proteogenic amino acids.

To elucidate the resistance of hybrid larch (Larix spp.) clones to Dendrolimus superans and their chemical defense mechanisms, 46 clones were subjected to field feeding-infestation experiments. Insect resistance was comprehensively evaluated by the membership function method from three perspectives: insect weight change, frass weight and defoliation rate (damage consequences); pre- and post-infestation physiological indices (constitutive and induced resistance). High- and low-resistance materials screened from each perspective were further analyzed by widely targeted LC-MS metabolomics. The three damage indices differed highly significantly among clones (P<0.01) and were positively correlated, indicating abundant genetic variation in resistance. Clones such as -30-C, -35-A and -21-A showed strong overall resistance, whereas -28-E and -20-A showed high constitutive and induced resistance. Defence enzyme activities (POD, SOD, CAT, PPO and PAL) correlated closely with secondary metabolites (total phenolics, tannins and total flavonoids) and nutrients, and rose synergistically after feeding, reflecting synergy between constitutive and induced resistance. In total, 3,577 metabolites were identified, predominantly terpenoids (>90%). The three groups yielded 1,010, 1,086 and 922 differentially accumulated metabolites, and PCA separated high- from low-resistance materials in all groups. These metabolites were significantly enriched in metabolic pathways, biosynthesis of secondary metabolites, monoterpenoid and diterpenoid biosynthesis, and phenylpropanoid biosynthesis. In conclusion, resistance of hybrid larch to D. superans is jointly conferred by the defence enzyme system, secondary metabolite accumulation and a 'terpenoid-phenolic' metabolic network: terpenoids and phenolics act synergistically via the terpenoid and phenylpropanoid pathways, together forming the core of chemical defense and supporting the breeding of insect-resistant clones.

Pectobacterium, a member of the soft rot Pectobacteriaceae (SRP), can cause highly destructive soft rot and blackleg diseases on many important crop and other host plants, particularly potato (Solanum tuberosum). Pectobacterium and its sister genus Dickeya have common virulence strategies, yet a highly diverse "open" pangenome containing many unique and hypothetical proteins. Here, we used randomly-barcoded transposon insertion-site sequencing (RB-TnSeq) on six Pectobacterium species to identify genes involved in bacterial survival and growth in potato tubers, as well as in minimal media conditions containing different single carbon sources. The strains evaluated varied broadly in genomic content as well as host plant and year of isolation. We focused on identifying conditionally essential genes across all strains tested, as well as genes that were important for growth in one or few of the tested strains. These results highlight the importance of many common genes involved in metabolism and maintaining cellular homeostasis, as well as shared and rare virulence traits including polygalacturonate catabolism, genes involved in polysaccharide biosynthesis, immunity proteins, and (for P. brasiliense) self-tolerance of the antibiotic carbapenem. This large dataset of conditionally essential genes, generated across Pectobacterium species during soft rot infection of the agriculturally important plant host potato, can serve as a rich dataset for further hypothesis generation and identification of important and/or novel in planta fitness factors.

Spontaneous vegetable fermentation is driven by the indigenous microbiota associated with the raw ingredients. The process often involves an initial heterofermentative phase followed by a predominantly homofermentative phase, both ensured by lactic acid bacteria (LAB). While microbial dynamics during fermentation are well characterized, phage communities remain poorly understood, despite evidence linking phage activity to LAB succession. We monitored viral and bacterial succession during sauerkraut and carrot fermentations using a metagenomics approach to characterize community dynamics. To further assess the impact of phage PLE001 on bacterial community dynamics, we developed a simplified carrot fermentation model comprising representative strains of Rahnella aquatilis, Leuconostoc mesenteroides and Lactiplantibacillus plantarum. Two consortia were assembled: one containing the PLE001-susceptible strain L. mesenteroides CIM3004, and the other containing also its phage-resistant derivative, CIM3004 BIM. We observed major shifts in the viral community during the first week of fermentation, coinciding with changes in bacterial structure during hetero- and homofermentative phases transition. In the simplified model, phage PLE001 reduced the population of the host during the heterofermentative phase but did not affect overall fermentation dynamics nor metabolite production. These results suggest that bacteriophages contribute to bacterial succession during fermentation, while communities remain resilient to phage infection, preserving fermentation stability.

Ensuring the safety and suitability of proteins considered for food and biotechnology applications is an important challenge. Computational approaches can facilitate the prioritisation of candidate proteins for further investigation, particularly given the rapid growth of genomic and metagenomic sequence datasets. In this study, we developed computational models to distinguish edible plant tissue-associated (EPTA) proteins from non-EPTA proteins. The positive class comprised 939 proteins with experimentally validated evidence of expression in plant tissues commonly consumed by humans, hereafter referred to as EPTA proteins. The negative class comprised 939 proteins annotated in UniProt as toxic, allergenic, or antinutritional, hereafter referred to as non-EPTA proteins. Models were evaluated using five-fold cross-validation with a cluster-based partitioning strategy, ensuring that no protein in the validation set shared more than 40% sequence identity with proteins in the corresponding training set. First, composition-based sequence features were evaluated to establish a baseline for EPTA or non-EPTA classification, with the Extra Trees (ET) classifier achieving a validation AUROC of 0.848. To incorporate evolutionary information beyond amino acid composition, PSSM profiles were subsequently evaluated, yielding an AUROC of 0.842 with an SVM classifier. Finally, to capture richer sequence patterns and contextual information using protein language modelling, we evaluated ESM2-based models. The fine-tuned ESM2-T33 model achieved the highest validation AUROC of 0.925 and was selected as the final predictive model. The resulting approach was implemented as the EdiProPred web server (https://webs.iiitd.edu.in/raghava/edipropred/) and standalone software for computational classification of plant proteins. The EdiProPred framework provides a useful computational tool for prioritising candidate plant proteins and facilitating the exploration of novel protein resources for food and biotechnology applications.